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Some of the arguments which support the strong concensus for an Q a = 0.3, A = 
0.7 model are reexamined. Corrections for Malmquist bias, local flow and mctallic- 
ity suggest a revised value for H of 63 ± 6 km/s/Mpc, improving the age problems 
for an Q = 1 universe. The latest CMB results may require a high baryon density 
and hence new physics, for example a strong lepton asymmetry. Difficulties for 
the fi D = 1 model with cluster evolution, the baryon content of clusters, and the 
evidence from Type la supernovae favouring low Q , A > models, are discussed 
critically. 



1 Introduction 

A strong concensus has developed in support of the cosmological parameter 
set: 

fl = 0.3, A = 0.7, H = 75, t a = 12.6 Gyr. 

Is this justified ? How certain is the evidence that A > ? In this paper I 
test how strong some of the evidence is, and show that the latest CMB results 
may point towards some interesting radical alternatives. 

Some of the problems for an f2o = 1> A = model are: 

• observed H t too high 

• Type la supernovae imply A > and CMB gives il + X = 1 

• cluster abundance evolution favours low o 

• baryon fraction in clusters favours low £l 

• P(k) for galaxies requires T ~ il a h ~ 0.2 

However these problems may not be as insuperable as they appear. 



2 Hubble constant, H a 

New distances from the HST Key program on Cepheids, and their application 
to Type la and II supernovae, and the Tully-Fisher, D n — a, and surface 
brightness fluctuation methods, have been reviewed by Freedman (2000) and 
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Mould et al (2000, and references therein). In these papers the approach is 
to determine the Hubble constant via different distance methods and then 
combine these estimates statistically. 

The philosophy of the present work is that of Rowan-Robinson (1985 
(RR85), 1988), to form weighted mean distances of groups and clusters and 
then restrict the analysis to high weight distances in order to minimise Malmquist 
bias 

An important ingredient in solving for the Hubble constant is the flow 
model. Estimation of peculiar velocities is one of the largest remaining uncer- 
tainties in H Q . Mould et al (2000) use a simple 3-attractor model. The PSCz 
flow model of Rowan- Robinson et al (2000), based on 88 groups, 854 clusters 
and 163 voids, gives much better resolution of the local flow. 

For Type la supernovae it is important to correct for internal extinction 
in the host galaxy. Phillips et al (1999) have carried out an important analysis 
of internal extinction for Type la supernovae and I have used these estimates 
where available, otherwise the RR85 prescription. I use the calibration of Gib- 
son et al (2000), based on 8 Type la supernovae with HST Cepheid distances, 
for which the absolute magnitude at maximum light is then < Mb > = -19.47. 

An important correction to the HST Cepheid distances is for the differences 
in metallicity between Cepheids in the LMC and the HST program galaxies. 
Mould et al (2000) estimate that this increases the HST Cepheid distance scale 
by 4%. 

The HST Key program is based on fi (LMC) = 18.5, but Feast (2000) 
estimates that this may be too low by 7%, based on Hipparcos parallaxes 
(after correction for metallicity effects in Cepheids). 

Figure 1 shows how the mean Hubble constant, calculated for groups and 
clusters with well-determined mean distances, varies with the assumed mini- 
mum weight. Inclusion of lower weight cluster and group distances leads to 
an increased mean Hubble constant, almost certainly the effect of Malmquist 
bias. The results are shown for the Mould et al (2000) 3-attractor flow model, 
and for the flow model derived from the PSCz analysis (Rowan-Robinson et al 
(2000). The difference between the flow models is not large, but it is significant. 

Figure 2 shows the Hubble diagram for groups and clusters with distances 
less than 400 Mpc. The mean Hubble constant for clusters with weight ( 
= TW/a 2 ) > 15 is found to be H Q = 63 ± 6 km/s/Mpc (Rowan- Robinson, 
2000a, in preparation), after correction for the effect of metallicity (but with 
Ho{LMC) = 18.5). I have adopted the uncertainty quoted by Mould et al 
(2000). 
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Figure 1: Effect of different minimum weights on mean value of Hubble constant. Filled 
circles: PSCz flow model (Rowan- Robinson et al 2000). Open circles: 3-attractor model 
(Mould et al 2000). No correction for metallicity has been applied. The error bars include 
only statistical uncertainties. 

2.1 Comparison with other estimates 

Cepheids: Frccdman (2000) have summarized the results of the HST Key 
Program using Cepheids in galaxies out to 20 Mpc to determine the Hubble 
constant and found H a = 74 ±8 km/s/Mpc. Mould et al (2000), with their 
3-attractor model for the local flow, find H Q = 71 ±6 km/s/Mpc. The dif- 
ferences between my estimate and Mould et al's can be assigned as follows: 
improved flow model AH a = -2, correction for Malmquist bias -2, correction 
for metallicity -3. 

Supernovae Type la: Branch (1998) has given an excellent review of 
the current situation. Initially widely disagreeing estimates of how the abso- 
lute magnitude at maximum light depends on decay rate have to some extent 
converged. Hamuy et al (1996) analyzed a large sample of nearby supernovae 
and find H Q = 63 ±4.5 km/s/Mpc. Saha et al (1997) analyzed 7 Type la su- 
pernovae with HST Cephcid distances and found H Q = 58 ±8 km/s/Mpc (61 
if they included a relationship between Mb and AM15). Branch et al (1996) 
found H a = 57 ±4 km/s/Mpc from a colour matched sample of supernovae and 
Tripp (1997) found H a = 60 ±5 km/s/Mpc from a sample matched according 
to A77J15. Riess et al (1996) used a template method in which multicolour pho- 
tometry is used to characterize supernovae in a 1-parameter sequence ('MLCS') 
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Figure 2: Hubble diagram for groups and clusters with weight > 3. 

to obtain H Q = 64 ±6 km/s/Mpc. 

Several groups have brought theoretical models to bear on the determina- 
tion of the Hubble constant using Type la supernovae. Hoflich and Khokhlov 
(1996) compared 26 supernovae with model light curves and found H = 67 
±9 km/s/Mpc. Branch (1998) suggests this should be revised to 56 ± 5. The 
same two authors found H Q — 55 if they included a theoretical version of the 
Mb - AMi 5 relation. Nugent et al (1995) fitted non-LTE model spectra to 
observations and found H a = 60 + 14,-11. In a recent analysis Tripp and 
Branch (1999) conclude that the best estimate for the Hubble constant from 
Type la supernovae was H = 62 ±5 km/s/Mpc. 

Gravitational lens time delay: 

An analysis by Falco et al (1997) of the gravitational lens time delay system 
0957+561 gave H Q = 62 ±7 km/s/Mpc. Koopmans and Fassnacht (1999) use 
5 gravitational lenses to determine H Q — 64 ±11 km/s/Mpc 

Sunyaev-Zeldovich effect: Recent work on SZ clusters includes Myers 
et al (1997), Birkinshaw (1999) and Reese et al (2000). Birkinshaw et al (2000) 
has given H = 54 ± 8 ± 10 km/s/Mpc for 9 clusters. 



To summarize, my estimate for H 0l 63 ± 6 km/s/Mpc agrees well with 
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independent estimates. 

With f2 Q = 1,A C = 0, the age of the universe becomes 10.2 ±1.0 Gyr, 
while for A = 0.7, o = 0.3, the age becomes 15.0 ±1.4 Gyr. These estimates 
can be compared with the estimate of Chaboyer et al (1998) for the age of the 
globular clusters, 11.5 ± 1.3 Gyr, to which must be added the time to form the 
globular clusters, 0.2-2 Gyr, depending on the redshift of formation. Neither 
cosmological model would appear to be decisively favoured (or ruled out) by 
the age estimate. 

3 CMB and n b h 2 

The pre-Boomerang and Maxima concensus was that f2 Q = 0.3, A Q = 0.7, Qbh 2 = 
0.019 ± 0.002, h = 0.68, n = 1. However this provides a very poor fit to the 
combined Boomerang and Maxima CMB data. In particular the second peak 
in the angular power spectrum is found to be too weak for this parameter 
combination. 

The concensus fix (Jaffc et al 2000) is to take a higher baryon density, 
fl b h 2 = 0.03. Unless there is a problem with Big-Bang nucleosynthesis calcu- 
lations, this would require the deuterium abundance to be significantly lower 
than even the low estimates of Buries and Tytler (1998a, b) (but see Tegmark 
et al 2000). It will be interesting to see whether the weakness of the second 
peak is confirmed in future observations. 

This approach also provides no physical explanation of why the vacuum 
energy density drops at the end of inflation by a factor 10 110 , only to achieve 
a significant dynamical role again at the present epoch. It is disappointing 
that quintessence models, which attempt to give some physical meaning to the 
evolution of A do not seem to fit the data as well as a constant A does. 

4 A radical alternative: lepton asymmetry 

A radical alternative is to consider a strong lepton asymmetry (Harvey and 
Kolb 1981, Kang and Steigman 1992, Lesgourgues and Pastor 1999, Kinny 
and Riotto 1999, White et al 2000, Lesgourgues and Peloso 2000, Esposito 
et al 2000). To maintain charge constancy, this has to be an asymmetry in 
the neutrino sector. The nucleosynthesis implications have been explored by 
Kang and Steigman (1992). The enhanced relativistic energy-density speeds 
up the expansion. Neutrons and protons decouple earlier, which enhances the 
abundance of neutrons, but this can be compensated by increased neutron 
decay due to an enhanced v e abundance. To get consistency with the ob- 
served primordial light element abundances, we then need a higher baryonic 
density. Lepton asymmetry also has the consequence (similar to the effects 
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of a positive A or a hot dark-matter component) of postponing the epoch of 
matter-radiation equality, thereby reducing the parameter T and improving 
the fit to the oberved P(k) relative to a pure CDM scanario. 

The lepton asymmetry can be characterized as an effective number of 
neutrinos, N e f / . The ratio of total energy density in extreme relativistic par- 
ticles to the energy density in photons, at the epoch of nucleosynthesis is then 
Per/ Pj = 11/4+ 7/8 N e ff. The expansion speed-up factor after e± anni- 
hilation can be written S a = [ 1 + 0.135 (N e ff - 3)] 1 / 2 . Lesgourgues and 
Peloso (2000) find a reasonable fit to the Boomerang and Maxima data with 
N eff = 6,n o = 1,A = 0,n b h 2 = 0.028, h = 0.70, n = 1 and r (reioniza- 
tion optical depth) = 0.05. Esposito et al (2000) find that the best fit to 
the Boomerang and Maxima data, retaining consistency with observed light 
element abundances, would require N v = 9 ± 4(2 — a). 

However to get a good fit to large-scale structure data, eg P(k), we need 

r = n h/s 2 = n h/[i + 0.135^ - 3)] ~ 0.2 

so for Q Q = 1, h = 0.63 we need N v ~ 19. 

The analysis of Esposito et al (2000) suggest that the best compromise 
would be N u = 13, which would give T = 0.26. This would give consistency at 
the 2 — a level with the CMB data, light element abundances and large-scale 
structure. 

A completely different radical alternative consistent with fl = 1,A = 
0, ilfc/i 2 = 0.019, and consistent with both the CMB data and large-scale struc- 
ture data is the proposal for a phase transition during inflation by Barriga et 
al (2000). 

5 The evidence for positive A from Type la supernovae 

There has been a great deal of excitement about the possibility that Type 
la supernovae yield positive evidence that A > 0. Recently two groups have 
published results on the Hubble diagram for Type la supernovae, with over 
100 supernovae now discovered at z > 0.3 (Schmidt et al 1998, Garnavich 
et al 1998, Riess et al 1998, Perlmutter et al 1999). Both groups claim that 
models with positive cosmological constant are preferred, and that models with 
A = 0.7, £1 = 0.3 provide the best fit to the data. The strength of the signal is 
that Type la supernovae at redshift 0.3-0.9 are about 0.25 magnitudes fainter 
than local supernoave, if an f2 Q = 1 Einstein-de Sitter universe is assumed. 
Claims that this is a 7-8 a effect therefore depend on a very precise homogeneity 
of Type la supernovae. The key element in reducing the scatter in Type la 
supernova absolute magnitudes at maximum light has been the correlation 
between absolute magnitude and decline rate (Ms — Amis), discussed above. If 
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Figure 3: Hubble diagram for SZ clusters (open circles, data from Myers et al 1997, Birkin- 
shaw 1999, Reese ct al 2000), Type la supernovae (filled circles, data from Pcrlmuttcr ct al 
1998) and gravitational lenses (crosses, data from Koopmans and Fassnicht 1999). 
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one looks at the paper by Hamuy et al (1996) where this relation is established 
for 29 local supernovae, one finds that the situation is not quite as impressive as 
has been presented. It would seem reasonable that to talk about a relationship 
between the absolute magnitude at maximum and the decline rate over the next 
15 days, it would be necessary to have detected the calibrating supernovae prior 
to maximum. In fact only 10 of the local supernovae were first observed at least 
one day before maximum. For these 10 there is indeed a Mb — Atoi 5 relation, 
but its significance is much reduced. If we derive the calibration from these 
10 local supernovae and apply it to the distant supernovae, the significance of 
the signal is reduced from the claimed 7-8 a to only 2-3 a, depending on which 
calibration is used (Rowan-Robinson 2000b, in preparation). It appears that 
the calibrating relation needs to be placed on a much stronger basis with nearby 
supernovae before it can be used to establish the reality of a cosmological 
constant. A good test of homogeneity would be to find several supernova in a 
high redshift cluster. 

Fig 3 shows the Hubble diagram, velocity versus luminosity distance, for 
snla, SZ clusters and gravitational lens time delay systems. The theoretical 
curves are shown for an assumed Hubble constant of 63 km/s/Mpc. The ev- 
idence from snla supporting A > comes from supernovae with distances > 
3000 Mpc (ie the fainter ones). The gravitational lenses, reaching to signifi- 
cantly greater redshifts, do not support A > 0. 

There are also some theoretical uncertainties. Since we do not know for 
certain whether nearby supernovae are due to white dwarf deflagration or to 
white dwarf mergers, there is the possibility that the proportion of these two 
types changes with epoch and this could affect the mean absolute magnitude. 
Hoflich et al (1998, 2000) also point out that uncertainties and evolution of the 
initial composition in supernovae can have a significant effect on the determi- 
nation of cosmological parameters using supernovae. 

6 Other problems: the cluster baryon fraction, cluster abundances 

• cluster evolution: strong negative evolution would be expected in the 
space density of rich clusters in an Q = 1 universe because growth of 
structure continues to the present day. This has been claimed not to 
have been seen by Bahcall et al (1997), Fan et al (1997), Carlberg et 
al (1997), Eke et al (1998), Bahcall et al (2000). Bahcall et al (2000) 
estimate that ft = 0.16 ± 0.05 and argue that fl Q = 1 is strongly ruled 
out. However using an X-ray selected sample of clusters, Blanchard et al 
(2000) estimate il = 0.74 ± 0.18, when account is taken of the cluster 
temperature distribution function. 
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• the cluster baryon fraction: Many groups have confirmed the original 
finding of White et al (1993) that the baryon content of clusters is too 
high for an f!„ = 1 universe ( White and Fabian (1995), Loewenstein and 
Mushotzky (1996), Mulchaey et al (1996), Evrard (1997)), if the baryon 
density is taken to be fl^h 2 = 0.019 ± 0.002. If however the higher 
baryonic density implied by the weak second CMB power spectrum peak 
seen by Boomerang and Maxima is correct, then the problem is not quite 
so acute. For example Birkinshaw (2000) estimates the fi&/0 o = 0.07, 
and if this is combined with ili,h 2 = 0.028, h = 0.63, we deduce that o 
= 1.008 ! 

In conclusion there are two viable scenarios consistent with the evidence 
on cosmological parameters: an Q — 1, A = universe or a Q — 0.3, A = 0.7 
universe. Although a number of lines of evidence, including high rcdshift Type 
la supernovae, favour the latter, the evidence from CMB fluctuations may 
require new physics to resolve. 
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